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Abstract 
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P±W- We find —dE/dz = \o. s N c p 2 LW holds independent of the details of the parton-nucleon 
scatterings so long as L is large. Numerical estimates suggest that p^-broadening and energy 
loss may be significantly enhanced in hot matter as compared to cold matter, thus making 
the study of such quantities a possible signal for quark-gluon plasma formation. 
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1 Introduction 



In a recent paper Q, hereinafter referred to as BDMPS II or simply as II, we investigated the 
problem of induced gluon radiation and radiative energy loss of a high energy quark or gluon 
traversing hot QCD matter of finite volume. This induced radiation may be an important tool 
for studying the physics of relativistic heavy ion collisions. In particular, contrasting hot and 
cold matter results may provide a signal for quark-gluon plasma formation. The induced gluon 
radiation may be measurable by studying the medium dependence of the characteristics of high- 
p±_ jet and hadron production. 

The present paper is devoted to the study of radiative energy loss and jet p_i_-broadening in 
nuclear matter. We shall show that these two phenomena are related in a model-independent 
way. Scatterings encountered by a fast parton in a nuclear medium may be described within a 
diffusion picture. The jet p^-broadening is determined by a classical diffusion equation. Energy 
loss is described in a picture where the fast parton propagation is treated classically whereas 
the induced gluon radiation is a quantum mechanical phenomenon. Interference effects play a 
crucial role in the production of gluons which determine the induced radiation spectrum. 

The Gyulassy-Wang model || used in II to describe multiple scatterings in hot matter is here 
replaced by the Glauber model || of multiple scattering in which parton-nucleus scattering is 
given in terms of independent scatterings off the nucleons making up the nucleus. 

The result which we find for radiative energy loss is similar to that for hot matter, namely 
that the total energy loss is proportional to the square of the length of the traversed nuclear 
matter, L. The coefficient of L 2 may be estimated from other, nonperturbative but measurable, 
quantities such as the small-x gluon distribution of the nucleon or, as discussed by Luo, Qiu 
and Sterman in Q, the transverse momentum which a jet, or dijet, receives as it passes through 
a nucleus. In this latter case, we find an interesting general relation between the radiative 
energy loss per length and the p^-broadening (defined as the characteristic width p 2 xw of the 
transverse momentum distribution) of a high energy parton passing through a nucleus, —dE/dz = 
^QZgNcP^yy. Numerical estimates suggest possibly large differences for pj_-broadening and energy 
loss for partons in hot matter as compared to nuclear matter. 

The outline of the paper is as follows. In section 2 the basic equations are established for 
jet p^-broadening and for the induced gluon radiation spectrum. These equations are solved in 
section 3 to logarithmic accuracy. Parameters appearing in the solutions are related to (other) 
phenomenological quantities, such as the gluon distribution of the nucleon, in section 4. Nume- 
rical estimates are given in section 5. In Appendix A, our basic equation for jet broadening is 
derived in detail. In Appendix B we relate the elementary parton nucleon scatterings, as they 
are used for p^-broadening and radiative energy loss, to the gluon distribution of the nucleon. 

2 The equations 

In this section we give the basic equations for jet transverse momentum broadening and induced 
gluon radiation for a high energy jet traversing cold nuclear matter. The equation for jet broad- 
ening is straightforward to write down and represents diffusion, in transverse momentum, due to 
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independent multiple scatterings of the jet as it encounters the various nucleons of the nucleus. 
The equation for induced energy loss has exactly the same form as found in II. In each case the 
fundamental assumption is that a parton is multiply scattered off uncorrelated scattering centres 
in the material. That is, we suppose that the Glauber picture of multiple scattering holds for a 
high energy parton passing through a nucleus. At extremely high energies the Glauber approx- 
imation breaks down due to coherent effects in inelastic particle production ||. We limit the 
energies which we consider in order to avoid this complication. It is an important goal for the 
future to better understand how severe these limitations are in practice. 



2.1 The basic elements of the equations 

We formulate the multiple scattering in close analogy with the procedure used in || within the 
Gyulassy-Wang model. Thus, we first need to find the cold matter quantity corresponding to 
the single scattering probability distribution. 

In jet p_i_-broadening the high energy parton receives a large transverse momentum due to 
the many scatterings in passing through a large nucleus. This means that in the individual 
parton-nucleon collisions all gluon couplings to the high energy parton occur with a small im- 
pact parameter separation between the parton in the amplitude and in the complex conjugate 
amplitude. This, perhaps surprising, result can be seen by referring below to ( |B.8| ) . If one 
inserts a complete set of scattering states between the two F^s in the right-hand side of that 
equation, the cross section is given explicitly in terms of a production amplitude times a complex 
conjugate amplitude. The amplitude F? 4 is evaluated at y = while in the complex conjugated 
amplitude y 2 = B 2 /fi 2 with fi a scale to be defined shortly. With such small impact parameter 
separations, B 2 oc L~ l , the one-gluon approximation becomes accurate. (The argument of the 
parton-gluon coupling is just the impact parameter separation.) 

In general, we allow breakup of the struck nucleon in the elementary parton-nucleon col- 
lision. The basic high energy quark-nucleon scattering amplitude is shown in Fig.l and the 
corresponding cross section is denoted by da/d 2 q. 




Figure 1: Basic quark-nucleon scattering amplitude. 
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A crucial assumption of the Glauber model we are using to describe multiple scattering is the 
fact that successive scatterings are independent. That is, the collision process at a given centre 
does not depend on the collisions which occur before or after. If the particles produced in the 
elementary collision in Fig.l have a formation time smaller than the average time between two 
scatterings (the mean free path A) then the independent scattering picture is clear. If that 
formation time is longer one expects corrections to independent scattering ||. In the following 
sections and Appendix B we shall relate an elementary scattering to the gluon distribution of 
the nucleon. 

An estimate of the size of the collective effects which undermine the Glauber model is the 
strength of shadowing corrections for the gluon distribution of a nucleus. Phenomenologically, 
these corrections are always less than a factor of 2 even for quite small values of x |7|||. Given 
this uncertainty, we expect the independent scattering picture to be a reasonable framework for 
studying jet p_i_-broadening and the induced gluon spectrum for high energy partons. 

Because da/cPq obeys diffractive kinematics, the four-momentum transfer squared, q 2 , is 
equal to q 2 the two-dimensional momentum transfer squared. It is then convenient to define 



as the total elementary cross section. The reader may worry about the meaning of a since the 
integration over small transverse momenta is not under perturbative control. We may imagine 
that the integral on the right-hand side of fl2.1|) has an infrared cutoff. The parton-nucleon cross 
section cannot be reliably calculated in perturbative QCD. However, we shall show that the 
cutoff dependence disappears (in the logarithmic approximation) when one calculates quantities 
insensitive to small momentum transfer scattering such as jet p^-broadening and induced gluon 
radiation. This cutoff independence is due to a cancellation of long-distance effects between the 
real inelastic production shown in Fig.l and forward elastic scattering, a virtual contribution. 
This cancellation is required by the probability conserving formalism which we use to describe 
multiple scattering. Any divergence in real production must be compensated by a corresponding 
virtual contribution (the no-scattering term on the right-hand side of ( |2.7|) ). This virtual contri- 
bution in our formalism shows up in the (cutoff dependent) mean free path for the quark given 



with p the nuclear matter density. 

In order to deal with dimensionless variables, and functions, it is convenient to introduce a 
(somewhat arbitrary) scale p 2 , with p representing a typical momentum transfer to the quark 
in a quark-nucleon collision. Then Q = q/p and U = p/p are the dimensionless transverse 
momentum transfer and transverse momentum, respectively. The probability distribution for 
the scattered quark is given by 





by 



A = [pa] 





(2.3) 



where we note that / d 2 QV(Q 2 ) = 1. The "potential" V, defined by (|2~3j ), will play exactly the 
same role in cold matter calculations as did a similar quantity, also called V, for hot matter |J. 
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As was done earlier for hot matter it is convenient to define the multiple scattering of a high 
energy parton in terms of the probability distribution for a single scattering, V(Q 2 ), and the 
probability distribution for distances A between successive scatterings [0,0, 

P(A)dA = e - A ' x< ^ . (2.4) 
A 

2.2 Equation for parton ^-broadening 

In calculating the transverse momentum given to a high energy parton by multiple scattering 
with the nucleons of a nucleus it is not necessary, in the leading approximation, to include 
the transverse momentum given to the parton by induced gluon emission since the emission 
is formally of higher order in a s . Thus the transverse momentum given to the jet in passing 
through a nucleus comes directly from the scattering of the projectile parton off the nucleons in 
the nucleus. 

Suppose a high-energy parton is produced in the nucleus in a hard collision at z = with an 
initial transverse momentum distribution fo(U 2 ) with 

d 2 Uf (U 2 ) = l. (2.5) 

(If the parton initially has no transverse momentum, then fo(U 2 ) = -5{U 2 ).) In terms of the 
scaled variable 

t={~ , \r = ^ (2.6) 
for a parton of colour representation R, the transverse momentum distribution at z — t\n is 

f(U 2 , t) = f {U*)e* + jf dtfe-W J d 2 Qf (0 - Q) 2 , f ) V{Q 2 ) . (2.7) 

The distribution f(U 2 , t) satisfies the normalization condition ( |2.5| ) for all t. We derive ( ^.7| ) using 
a multiple scattering formalism in Appendix A. The first term on the right-hand side of (|2.7|) 
corresponds to no interaction with the nuclear medium where the e~* represents the probability 
that no interactions have occurred in a distance z = A#£. The second term on the right-hand 
side of ([2.7D corresponds to a final scattering at t' before the parton reaches t with transverse 
momentum U. One can recast ( |2.7|) as a differential equation by taking a "time" derivative on 
both sides of that equation. One finds 

df(U 2 ,t) £fTT2 . , . ^ 2 ^ rl ^ 



, n -f{U\ t) + J d 2 Qf((U - QY, t)V(Q 2 ) . (2.8) 

A different, but equivalent, form of ( |2.8| ) for /, which makes its interpretation as a kinetic master 
equation more transparent, is 

8f{U J t lt) = ~ J f(U\ t)V((U - U') 2 )d 2 U' + I f(U' 2 , t)V((U' - Uf)d 2 U' . (2.9) 



5 



When partons pass through the distance dt the first term on the right-hand side (loss term) 
accounts for partons which are scattered out of the direction [/,[/—> U', weighted by the 
(normalized) scattering cross section V. The second term (gain term) counts those partons which 
are scattered into the direction U from all other directions U', U' — > U. One can diagonalize 
( p78|) by defining 

f(B 2 ,t) = J d 2 Ue- iS - ff f(U 2 ,t) , (2.10a) 

and 



V(B 2 ) = / d 2 Qe~ iS4 V{Q 2 ) , (2.10b) 



in which case 



where 



df(B 2 ,t) l D2r . /D2w ~ /D2 



dt 4 



4 



-B 2 v(B 2 )f(B 2 ,t), (2.11) 



v(B 2 ) = —(l-V(B 2 )). (2.12) 

Note that V^(0) = 1 because of our normalization of V(Q 2 ) as a probability distribution. 

Instead of searching for an exact solution one often approximates (|2.8|), in case of soft "po- 
tentials", by expanding f((U — Q) 2 , t) inside the integral through second order in Q. This leads 
to the diffusion equation [pf 



df(U 2 ,t) 1 



()f A q uV tjf(U 2 ,t), (2.13) 

where the "transport coefficient" q v is defined by 

% = Jd 2 QQ 2 V(Q 2 ). (2.14) 

The transport coefficient with respect to the longitudinal variable z (see ( |2.6|) ) and the transverse 
momentum q is then given by 

q = ^q v . (2.15) 



In cases where the integral in Q2.14D exists one finds that 

% = m- (2.i6) 



For example, if V were gaussian, V(Q 2 ) = Q , (|2.7|) and ( |2.14|) would be completely equivalent 



to standard treatments of parton transverse momentum broadening in nuclei [|T(| with /i 2 the 
average transverse momentum in an elementary collision. 



6 



However, in QCD ( |2.14| ) diverges logarithmically since V(Q 2 ) oc 1/Q A at large Q 2 . Corre- 
spondingly, v(B 2 ) has no finite limit for B 2 — > 0. Nevertheless, v(B 2 ) may be related to the first 
moment of V(Q 2 ) in a logarithmic approximation. To see this we use ( |2.10b| ) and ( |2.12| ) to write 



v{B 2 ) = — / d 2 Q 1 - e 



Doing the angular integral gives 

v(B 2 ) 



B 2 



47r r°° 



-iB-Q 



V{Q 2 ) . 



B 2 Jo 



dQ 2 {\- J (BQ))V(Q 2 



(2.17) 



(2.18) 



This integral is convergent in the ultraviolet. For BQ <C 1, we have 1 — Jq(BQ) jB 2 Q 2 so 
that there is a logarithmic region of integration in ( 2.18Q extending up to Q 2 ~ 1/B 2 . Thus, in 
the logarithmic approximation, for small B 2 



v(B 2 



1/B 2 



dQ 2 V{Q 2 )Q 2 



(2.19) 



2.3 Equation for the induced energy spectrum 

Because of the close correspondence between the multiple scattering in hot matter, as discussed 
in II, and the present circumstance we need not rederive the equation for the induced radiation 
spectrum. We may simply take over equations (4.16), (4.24) and (5.1) of that paper which read 







dr 

_ 3oz s Cr 
dujdz 27r 2 A„ 



f(U, r) 



1 - iW 2 )f{U, r) + / d 2 Qf(U - Q, t)V(Q 2 ) 



udl 



2Re li 



TO 



dr ( 1 

to 



r\ rd 2 B 1-V(B 2 



2tt 



5 2 



-B-f(B, 



k=0 



and 



SV | - -B 2 v(B 2 ) 



VK 



f(B,r); f(B,0) = --Bv(B 2 ) 



respectively. In (ggg)-(ggg) 

t = z/X, r = L/X 



2C F 



k , A 



2C Fx _ o A 



(2.20) 
(2.21) 

(2.22) 
(2.23) 



Equations ( |2.21| ) and ( p.22| ) give the induced gluon radiative energy spectrum for high energy 
partons having colour representation R. The importance of the parameter k = Xfi 2 /2uj has been 
discussed in 0. 



3 Energy loss and jet broadening 



In this section, we give the solutions to fl2~20D - flpp . We then compare the resulting formulas 
for dE/dz and p\w and note a simple relationship between these two quantities. 



7 



3.1 Jet transverse momentum broadening 

as 



It is straightforward to solve ( 2.11 



f(B 2 ,t) = f(B 2 ,0)e W {~B 2 v(B 2 )t} . (3.1) 



In this paper we only consider media whose length L obeys L/X ^> 1. Thus we need only 
consider t large in which case only small values of B 2 give a result which is not exponentially 
small. If v(0) were finite we could choose fi 2 to be the mean momentum transfer squared per 
collision so that v(0) = 1. In this case the transverse momentum distribution, f(U 2 ,t), would 
be a gaussian in U with mean transverse momentum squared (pj_) = fi 2 t. However, we expect 
a logarithmic dependence of v(B 2 ) for small B 2 because of the approximate scale invariance of 
QCD. We suppose that B 2 -^ \nv(B 2 ) <C 1 for small B 2 in which case it is clear from (|3.1|) that 
the values of B 2 for which f(B 2 ,t) is not small are 

B 2 < — i— . (3.2) 

Using /(0, 0) = 1, coming from the normalization of f(U 2 ,t) as a probability distribution in 
transverse momenta, we get 

f(B 2 ,t)*exp{-±B 2 tv(l/t)} . (3.3) 
Transforming back to [/-space gives 

?' a, ^~^M-T§m}- (3 ' 4) 



In neglecting the B ' = singularity of v in ( |3.4| ) we have lost the large- £/ tail of f(U , t) which 
is 

f{U\t) ~ V{U 2 )t ~ ±. 

The gaussian approximation ( p.4| ) is valid up to [7 2 -values where the exponential is as small as 
1/t <ti 1. This means that to evaluate the characteristic width of the distribution one can use 
). Defining p 2 LW as the value at which the distribution falls to 1/e of its peak value, we have 



pl w = V 2 J d 2 UU 2 f^\U\L/X R ) , (3.5) 

which results in 

Plw = fv(X R /L) L . (3.6) 

The linear growth of typical transverse momenta squared with L is expected and has been used 
for some time to discuss p_i_-broadening of high energy partons in nuclei |T0|,P 
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The individual parts on the right-hand side of ( |3.6| ) which multiply L are not uniquely defined. 
Nevertheless, we note that the quantity ^-v(Xr/L) is independent of the scale /i 2 and of any 
infrared cutoff we may have used to define a. Indeed v scales as 1/yU 2 so that [i 2 v is independent 
of /i 2 . More importantly, the factors of a cancel in v/X, and the dangerous small-Q region of 
integration is suppressed for small B in taking the difference indicated in ( |2.17| ) leaving the cutoff 
only as a scale for logarithmic dependence. Indeed, from ( |2.19| ) 



A W dQ W Q - (37) 

which explicitly shows the suppression of long-distance contributions. 

It is worth mentioning that the basic evolution equations ( p.llj ) and (|2.22| ), and thus their 



solutions, are also cutoff independent for large z if we substitute t = z/Xr everywhere in the 
equations and view / as a function of the physical variable z rather than t. We prefer to 
continue working with the variable t, however, as the correspondence with the procedure used 
for hot matter is more transparent in this variable. 



3.2 Energy loss 



Since ( |2.21 ) and ( |2.22| ) are equations identical to those encountered in II we may immediately 
write down the solution to ( 2.22|) as 



f(B,r) 



2mB d 



and we may write ( p. 21 



as 



B 2 dr 

udl 6a s Cji 
dudz irL 



exp •{ — —m ujqB tanc^oT } , 



In 



sin lu t 



where 



m 



1 

Ik 



ikv . 



(3i 



(3.9) 



(3.10) 



The argument of v is B 2 in (|3.8| ) and yk/\n(l/k) in ( p.9|) . Exactly as in II one may give an 
analytic expression for dE/dz, the energy loss per unit length, as 



dz 



and the total energy loss as 



dE as ° R ^ v(X/L) L = ^^ v(X/L)L 



X, 



x f 



AE = ^ £v(\/L) L 2 
8 X g 



(3.11) 



(3.12) 



In the above, a s = a s (k 2 ) is evaluated at a large momentum scale k 2 proportional to L. Indeed, 
as shown in II, the effective value of the impact parameter B 2 is of order X/L. 
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3.3 Relationship between energy loss and ^-broadening 



Comparing (|3lf ) and ( |3.11[) we note that 

dE a s N c 2 

- = ~r p±w (3 - 13) 

relating the energy loss per unit length with the typical transverse momentum squared that a jet 
receives in passing through a length L of nuclear matter. It is interesting that the relationship 
between these two quantities is completely independent of the dynamics of the multiple scattering, 
at least in the framework of our approach to these problems. Thus (|3.13|) holds equally well 
in finite length hot matter as well as cold matter. Perhaps even more surprising the coefficient 
relating p\ w to dE/dz is independent of the nature of the high energy parton passing through the 
hot or cold matter. Equation ( |3.13|) makes more precise the bound —dE/dz < \p\ w suggested 
sometime ago by Brodsky and Hoyer WL\ on the basis of the uncertainty relation. 



4 Determining tj-v, estimates 

2 

In this section, we shall discuss how can be determined phenomenologically and we shall give 
a crude estimate of this quantity in terms of the nucleon gluon distribution. 



4.1 Determining phenomenologically 

The broadening of the transverse momentum spectrum of high energy partons in nuclei has been 
widely discussed in the literature. In /z-pair and J/V'-production in high energy proton-nucleus 
collisions it is the incoming high energy quark (or gluon) which multiply scatters in the nucleus 
and in so doing broadens the /i-pair and J/ ip transverse momentum spectra as compared to the 
corresponding production in nucleon-nucleon collisions. Transverse momentum broadening has 
also been studied for outgoing jets produced in high energy proton-nucleus and photon-nucleus 
collisions. 

The phenomeno logical situation is somewhat confused at the moment. It would appear that 
outgoing partons receive much more transverse momentum due to multiple scattering in a nucleus 
|12|| than do incoming partons [13]. Theoretically, this difference is not understood PJT^j. It is 



not the purpose of this section to give a comprehensive study of the phenomenology of multiple 

2 

scattering of high energy partons in nuclei but rather to show how can be determined, so we 
limit ourselves here to relating our approach to that of Luo, Qiu and Sterman [|J who studied 
the p_i_-broadening of dijets produced in photon-nucleus collisions at Fermilab. 

Luo, Qiu and Sterman consider a single rescattering in the nucleus of dijets produced in high 
energy photon-nucleus collisions. They find that the transverse momentum squared given to the 
dijet by summing over all single rescatterings in the nucleus is (see equation (21) in |4]]): 

kUR)) = ^ 2 a s (Q 2 ) A 1 ? 3 \l QS . (4.1) 
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Here Q is the relative transverse momentum of the two jets making up the dijet, A is the atomic 
number of the nucleus and X 2 L q S is a dimensional parameter which characterizes the momentum 
transfer squared in a single collision. An expression for \\q S is given in J3J in terms of a new 
QCD nuclear matrix element, but the actual determination of A|q 5 is made by comparing (|4.1|) 
to experiment. A value \\qs ~ 0.05 — O.lGeV 2 is found. 
We may compare (|3.6| ) to ( j4.1|) by noting the following: 

(i) In the Luo, Qiu and Sterman calculation, the dijet consists of a high transverse momentum 
quark and a gluon with the quark and gluon transverse momenta nearly balancing. It is the 
nuclear contribution to the imbalance, namely, the transverse momentum of the whole dijet, 
which is expressed in (|4.1|) . Since the quark and gluon have very high transverse momentum 
they form a (colour triplet) system which is compact in transverse coordinate space. Thus, in 
a relatively soft rescattering with a nucleon in the nucleus this system should act like a single 
quark and this indeed is the case in the calculation done in 01] . We may then expect that our 
calculation of jet broadening of a single parton should apply to the dijet system. 

(ii) Luo, Qiu and Sterman only consider a single rescattering while we have summed arbitrary 
numbers of rescatterings in arriving at ( |3.6|) from (|2.7| ). However, in a multiple scattering, where 
the scatterers are uncorrelated, the resulting transverse momentum squared is just the sum of 
the transverse momenta squared of the individual scatterings. In 0J only a single scattering is 
explicitly done but (|4.1|) results from summing over all possible single scatterings. Thus (|3.6|) 
and (|4.1|) should be directly comparable. However, here we do sidestep the difficult issue as to 



exactly which (/c|,)-average has been done experimentally in evaluating X L q S and how close that 
average is to our p± w . 

The last piece of information we need is to relate L to the nuclear radius R. Since the dijet is 
produced in a hard collision it is produced uniformly in the volume of the nucleus. After being 
produced the average length of material that the dijet passes through is L = With this 
identification of L we may equate (|3.6|) to fl4.1|) and obtain 



/^ 2 , M/n 16 2 ^ 2N ^ 1/3 x2 



^v(X/L) = -^a s m — X L q S ■ (4.2) 



We shall give a numerical estimate of from (|4.2|) a little later on. 
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4.2 Relating ^yv to the gluon distribution 



One can get a crude relation between and the nucleon gluon distribution [15] and that is 
our purpose in this section. Indeed, Luo, Qiu and Sterman relate their parameter X\q S to the 
average colour field strength squared that a parton sees as it passes through a nucleon in the 
nucleus. In order to relate to the gluon distribution we use ( |2.19|) which gives 



v(B 2 



1 



dQ 2 Q 2 ^ 

n dQ 2 



m 2 /b 2 da R 

dq q -ty 

o dq 1 



(4.3) 



11 



where now da R /dQ 2 refers to the scattering of a parton of colour representation R with a nucleon. 
But in the one-gluon exchange approximation 



^ /B2 . 2 2 d(J R 

u dq z 



(We postpone until later on the prickly question of what value of x we intend to use in 



(4.4) 



We 



do note however that in the region that (4.4) will be used, a region of small but not too small 
x, the quantity xG(x, /i 2 / B 2 ) should have little dependence on x.) The simplest way to see that 
Q4.4D should be true is to consider the case where the target and jet are both quarks. In this case 



to lowest order in a, while 



da 
dq 2 

xG q (x,fi 2 /B 2 



2%C F a 



2\2 



N c (q 



QJeCt 



7T 



1 f 

ln T 2 



(4.5) 



(4.6) 



Using (|4.5|) in the left-hand side of (|4.4j) gives the right-hand side of that equation after using 
([4.6[), thus fixing the 2n 2 a s /N c factor on the right-hand side of (fOj), for a quark jet, and a factor 
4:TT 2 a s Cii/(N 2 — 1) for a jet of colour representation R. (A more complete, and perhaps more 
satisfactory, derivation of (\i.4\) is given in Appendix B.) The a s in (|4.4j) should be evaluated at 
the scale n 2 /B 2 as B/fi is the short distance over which the scattering occurs. Using ( [4.4[ ) in 



1|) along with = [pu R ] gives 



^v(B 2 ) 



A 



N 2 - 1 



pxG(x, p 2 /B^ 



(4.7) 



5 Some numerical estimates 



The basic results for p±— broadening and energy loss are given in (|3.6| ) and in ( |3.11 ), ( 3.12 ) 

2 

respectively. The parameter controlling these quantities is v, for which we have found the 



expressions ( |4.2[) and (|4.7| ). Our task here is to make rough estimates of this quantity. 

(i) In II, for hot matter having T = 250 MeV, we took /i 2 / A = 1 GeV/ fm 2 and, for L 
v w 2.5 giving 



10 fm, 



We get 



0.5 GeV 2 / fm. 



5 GeV 2 



from (|3.6|) and, taking a s = 1/3, 



-AE ~ 30 GeV 



10 fm 



L 
10 fm 
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for a quark jet. These are rather big numbers and one should not take their exact values too 
seriously. However, they do suggest that hot matter may be rather effective in stimulating 
radiative energy loss and in broadening the p_i_-distribution of high energy partons. 



(ii) For cold matter from (|4.7| ), with Cr = Cp = 4/3, using p = 0.15 fm one finds 

^ « — • a s ■ \xG(x)} GeV 2 /fm . 
A 25 L V n ' 

Taking a s = 1/2 and xG = 1 results in much smaller values 

L 



■ 2 ~ 0.2 GeV 2 



P±w 



10 fm 

and 

—AE » 2 GeV 



L 
10 fm 

These numbers may be reasonable. In particular, dp]_ w /dz = ^ ~ GeV 2 / fm is in agreement 
with the pj_-broadening of the /i-pair spectrum[13] so —dE/dz w 0.2 GeV/ fm(L/10 fm) may be 
a sensible, if small, estimate of radiative energy loss for high energy quarks. 

(iii) From (fO|), with X 2 LQS = 0.05 GeV 2 , one finds ^ ~ 0.8a s (Q 2 ) GeV 2 / fm. Taking 
a s (Q 2 ) = 1/3 gives 



A 

and leads to 

2 



2 ~ 

^ V ~: 0.3 GeV 2 / fm 



pl^^SGeV^ 



10 fm 

and 

— A_E ~ 15 GeV 



L 
10 fm 

These numbers seem too large. (They would correspond to xG ^ 20 — 25 if ( |4.7[ ) is used.) The 
problem here, as alluded to earlier, is that present experiments find a large p^-broadening and 
energy loss for outgoing partons giving A|q 5 ~ 0.05 — 0.1 GeV 2 while much smaller numbers are 
found for the broadening of the /z-pair spectrum. 
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A 



Appendix A. Basic equation for jet broadening 

Here we give a brief derivation of ( |2.7|) . The probability density, p n (U 2 , t), for the incoming parton 
to have a transverse momentum squared U 2 at time t after n collisions at times t X) t 2 , ■ ■ ■ t n and 
involving fixed transverse momentum transfers Qi, Q2, ■ ■ ■ Q n is f| 

(n \ n 

0-T.Qi) 2 II [viQ 2 ^-^] e-(*-*») . (A.l) 
i=i / e=i 

Thus, 

Cn \ n 

(u-j2Qi) 2 mm 2 ), (a.2) 
i=i / e=i 

where t = 0. 

The probability density, P n (U 2 ,t), to have momentum U at time t is obtained by integrating 
over the Qi and ti leading to 

P n (U 2 , t) = / n dt» f " dVi • ■ • r dt i Pn(U 2 , t) . (A.3) 

J fj[ Jo Jo Jo 

The momentum distribution / is given by the sum over all possible number of interactions, 

00 00 
f(U 2 , t) = £ P n (U 2 , t) = P (U 2 , f) + £ P„,(f/ 2 , t) , (A.4) 

n=0 n=l 

where the probability to have no interactions with the medium between times and t is 

Po(U 2 ,t) = f (U 2 )e- t . (A.5) 
It is easy to relate P n to P n -\. Inserting ( |A.2| ) into ( |A.3[ ) gives 



n-l 



P n (U 2 , t) = f d 2 Q n V(Q 2 n ) f dt n e-^ J] d 2 Q e V(Q 2 ) f " cft^ • ■ • / ' 

J Jo £^ JO JO 

• e-*"/of(^-4-E^) 2 V ( A - 6 ) 



=1 



Relabeling Q n as Q and t n as t', 

P n (U 2 ,t) = [dt'e-W J d 2 QV(Q 2 )P n ^((U~Q) 2 ,t') , (A.7) 
which has an obvious probabilistic interpretation. Finally, using (A.7) in (A.4) leads to (|2.7| ). 
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B 

Appendix B. Relating elementary scatterings to the gluon distribution 



In this appendix we derive (|4.4| ). Much of our discussion here follows [|J. We assume a one- 
gluon exchange approximation at the outset (see Fig.l), the justification of that assumption being 
the small coupling we find for the coupling of the exchanged gluon to the incident high energy 
parton. In order to have expressions for the gluon distribution which conform to conventional 
notation we suppose that the incident parton is moving in the negative z-direction so that p = 
and, using light-cone variables, p_ ^> p + . (In order to fix (approximately) the value of x in 
( |4.4| ) we must consider the parton propagating through the medium to be virtual though, in our 
covariant gauge discussion which follows, we may continue to assume that the exchanged gluon 
couples in an eikonal manner.) We suppose the nucleon, labeled by P, is at rest. Then 



da 
d 2 q 



<xg 2 / dq+dq-d?y(P\AZL(y)Ai.(0)\P)6((p + q) 



(b.i; 



with the gluon field with colour index a. The process is illustrated in Fig. 2 where and y are 
the coordinates where the gluon attaches to the high energy parton in the production amplitude 
and complex conjugate amplitude, respectively. 



p 


y p 


i 


1 — t~ 

°i 

i q 




p 


p 



Figure 2: Diagrammatic representation of the cross section (B.l). 
With (q + p) 2 ~ q 2 +p 2 + 2p_q + the q + — integration is easily done using the 6— function in (B.l): 



da 
d 2 q 



oc g 2 dq_d 4 y(P\A a Jy)A a jO)\P)e 



-iq+y--iq~y++iq-y 



(B.2) 



with q + m —{p 2 + q 2 )/2p^ m —p 2 /2p_ in ( [B.2|) . (We shall see later that q 2 <p 2 -) The q. 
integration gives a 5-function in y + (because of our large-p_ limit) so we find 



da 
d 2 q 



cx g 2 / d 2 ydy-(P\A a + (y + = 0,y^,y)A a + (Q)\P)e-^ y - +i ™ 



(B.3) 
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Now multiply both sides of ( |B.3| ) by q 2 . The q 2 on the right-hand side of ( |B.3| ) can be written 
in terms of derivatives of the A's giving 

f^9 2 j d 2 ydy4P\ ±A a + (y + = 0,y^y)^-A a + (0) \P) e"^^. (B.4) 

At lowest order in g 

F? + (y) = ^Al(y) - . (B.5) 

The second term on the right-hand side of (|B.5| ) is small when used in (|B.4j ) because the ?/_- 
derivative can be integrated by parts to give a q + while the + index on A + in the first term 
naturally yields a factor of P + when used in QB.4Q . (We shall shortly identify x = q+/P+ ~ 
—p 2 /2p_P + = —p 2 /s <C 1 with the Bjorken x- variable of the gluon distribution.) Thus, 

oc g 2 J d 2 ydy.(P\ F^(y + = 0, y_, y)F^(0) |p) e -*»^-W. ( B .6) 

Integrating flB.6j ) over q 2 between and fj, 2 /B 2 gives 

^ 2 /s 2 de- 
ep g 



J o d 2 qq 2 —^g 2 J Q dpJ^p) J dy_ e - teP +»-(P|J^ 4 (y + = 0, y_, $ -i^(0)|P) 



(B.7) 



where we have used zJo(z) — d/ dz{zJ\(z)). Values of p on the order of 1 dominate the integral 

(B.8) 



in ( B.7 ). Therefore in the leading logarithmic approximation we may write 



Jo ^^ik^ 9 ] ^ +y ~( P \ F Uy+ = 0,y-,y)FU0)\P) 



^2 =B 2/ fJ _2 



The integral on the right-hand side of ( |B.8| ) is proportional to a standard expression for the 
gluon distribution xG(x, fi 2 /B 2 ) at a momentum scale /i 2 / 'B 2 JT|]. The g 2 in ( |B.8|) (the a s in 
( |4.4D ) should also be evaluated at a scale fi 2 / ' B 2 since this is the limit on the q 2 integration set 
by the hard scattering. Once we know that the left-hand side of ( |B.8| ) is proportional to xG, the 
constant of proportionality is easily set using (|4.5|) and ( f4.6|) . 

Finally, we come to the value of x to be used. To determine the approximate value of x 
we must decide what value of p 2 should be used for the incident parton. We suppose that the 
incident parton (jet) is produced in the nucleus in a hard collision with p = 0,p- fixed and 
p + varying with the distance from the production point as p + ~ 1/z. (We recall that due to 
the uncertainty principle one cannot specify p + more accurately in a distance z.) Then for our 
elementary scattering 

^IP 2 !^, (B.9) 
where L is the length of the nuclear material and A# the mean free path. Since 

f<^~ 2 ~ -^v{\ R /L)p? (B.10) 



R 
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(see(3.2)), 



W\ ~ ^ (Rn: 



a quantity which we assume not to be large. (Note that for L = L cr , with L cr defined in jl]J 
one has /i 2 L 2 /A#p_ ~ 1.) Thus, we are justified in dropping q 2 compared to p 2 in determining x 
from (q + p) 2 = 0. 

From ( p.9|) we come to the conclusion that x should be chosen in the range 

m ix ^ik' (B12) 

with M the nucleon mass. 
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